Noname manuscript No. 

(will be inserted by the editor) 



Against all odds? Forming the planet of the HD 196885 
binary 

p. Thebault 



Received: date / Accepted: date 



Abstract HD196885Ab is the most "extreme" planet-in-a-binary discovered to date, 
whose orbit places it at the limit for orbital stability. The presence of a planet in 
such a highly perturbed region poses a clear challenge to planet-formation scenar- 
ios. We investigate this issue by focusing on the planet-formation stage that is ar- 
guably the most sensitive to binary perturbations: the mutual accretion of kilometre- 
sized planetesimals. To this effect we numerically estimate the impact velocities dv 
amongst a population of circumprimary planetesimals. We find that most of the cir- 
cumprimary disc is strongly hostile to planetesimal accretion, especially the region 
around 2.6AU (the planet's location) where binary perturbations induce planetesimal- 
shattering dv of more than Ikm.s Possible solutions to the paradox of having a 
planet in such accretion-hostile regions are 1) that initial planetesimals were very 
big, at least 250km, 2) that the binary had an initial orbit at least twice the present 
one, and was later compacted due to early stellar encounters, 3) that planetesimals did 
not grow by mutual impacts but by sweeping of dust (the "snowball" growth mode 
identified by Xie et al., 2010b), or 4) that HD196885Ab was formed not by core- 
accretion but by the concurent disc instability mechanism. All of these 4 scenarios 
remain however highly conjectural. 

Keywords Planetary systems ■ Binary Stars 



1 Introduction 



1 . 1 Planets in binaries 



Studying planet formation in binaries is of fundamental ir nportance, as a majority 
of ma in sequence stars are members of multiple systems (IDuquennov and Mavon 



199 Ih . Moreover, planets in binaries are no longer theoretical concepts, as about 
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20% of all known exoplanets have been found to iiihabi t multiple stellar systems 
( Desidera & Barbieril 2007 : Mugrauer & Neuhauserl 2009 1. Most of the planet-bearing 
binaries have large separations, often in excess of 1000 AU, for which the influence of 
the companion star on the planet region, and the planet formation process, is probably 
limited. However, a handful of planets inhabit much tighter binaries, which separa- 
tions as small as ~ 20 AU: G186, with a pl anet at 0. 1 1 AU and a companion at 18.4 AU 
JOueloz et al.1 12OOOI: Eagrange et al.1 Eood) . HD 41004 with a planet at 1.64 AU and 
companion at 23 AU dZucker et al.[ '2004*) a nd y Cephei, for which the plane t lies a t 



2.04 AU and the companion at 20.18 AU (Hat zes et aUl2003l:lNeuhauser et al.Ll2007h . 



For these systems (especially 41004 and j Cephei), the closeness of the companion 
must have had an influence on the way the planet formed and dynamically evolved. 
Studying how such planets came about is of great interest, all the more because close- 
binaries can be used as a test bench for planet formation models, by confronting them 
to an unsual environment where some crucial parameters might be pushed to extreme 
values. 

Historically, the first issue that has been investigated is that of the long term sta- 
bility of planet ary orbits in binaries. T he reference work on this issue remains proba- 
bly that of Hol man & Wiegert (Il999l) . who derived empirical expressions for orbital 
stability as a function of binary semi-major axis a^, eccentricity and mass ratio /i. 
Later studies have show n that, reass uringly, all known exoplanets in mu ltiple systems 
are on stable orbits (e.g. Dvorak et a l.. 2003; Haghighi pour et al. . 2010l) . although the 
case for HD41004 is not full y settled yet, as it depends on the yet unconstrained ec- 
centricity of the binary orbit (IHaghighipour et al.i I2OIOI) . 



1.2 Planet formation in binaries 

The next step is to study under which condition such planets can form, as the con- 
straints for planet formation might be very different from those for orbital stability. 
This issue is a much more difficult one. Planet formation is indeed a very complex 
process, believed to be the succession of several stages (e.g. iLissauer [ 1 19931) . each 
of which could be affected in very different ways by the perturbations of a secondary 
star. Not surprisingly, the effect of binarity on each of these different stages is usually 
investigated in separate studies. 

The initial phase of planet formation, i.e., the formation and evolution of a gaseous 
prot oplanetary disc , has b een investigated early on by Artvmowicz & Lubow (1994 
and ISavoniie et al.l (1 1994 ), who have shown that the circumprimary disc is tidally 



truncated by the companion. This truncation occurs at a location comparable to the 
outer limit for dynamical stability and thus safely lies beyond the position of all de- 
tected exoplanets. It could nevertheless pose a problem, especially for giant planet 
formation, as it deprives t he disc from a larg e fraction of its mass. For the specific 



case of y Cephei, however. iJang-Condell et a l. (2008) have found that there is proba- 



bly enough mass left in the truncated disc to form the observed giant planet. But an- 
other, potentially more troublesome consequence of disc truncation is that it shortens 
the viscous lifetime of the disc, and thus the timespan for gaseous planet formation. 
This effect seems to have been observationally confirmed bv lcieza et al .I (l2009h . who 
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found that young binaries with separation < 100 AU have a lo wer probably of h osting 
circumstellar dust in the innermost few AU around each star touchenellioiol) . even 
if some close binaries do show signs of a hot circumprimary disc. 

The next stage of planet formation, the condensation of small grains and their 
growth into larger pebbles and eventually kilometre-sized planetesimals, has not been 
extensively studied in the context of binary systems. One main reason is probably that 
this stage is the one that is currently the lea st understood even in the "normal" context 
of single stars (e.g. Blum & Wurmi 2008 ). so that extrapolating it to perturbed bina- 
ries might seem premature. A noteworthy exception is the study by Nelson (200d|) 
showing that for an equal-mass binary of separation 50 AU, temperatures in the disc 
might stay too high to allow grains to cendense. But, as ackowledged by the author 
himself, these results are sti ll prelimina r y and a full study of this issue has yet to 
be undergone. Very recentlv. IZsom et aL ( 201 1 ) showed that, even if grains can con- 
dense, binary perturbations might impend their growth by mutual sticking because of 
too high impact velocities. 

The stage for which the influence of a companion is probably best understood 
is the final step of planetary accretion, leading from Lunar- sized embryos to fully 
formed planets. Several studies have shown that the regions wh ere embryo accretion 



can proceed roughly correspond to those for orbital stabihty (IBarbieri et al.l 12002 



lOuintana etaP. l2007t lOuedes et all 120081 lHaghighit?oure t aLllioiO l). This is a fur 
ther reassuring result for all known exoplanets-in-bina ries, in particula r the ar chetypal 
case that is y Cephei, for which the extensive study of lKlev & NelsonI ( |2008|) showed 
that, // embryos can form around 2 AU from the primary, then they can evolve to 
form planetary cores at the present planet location. 



1.3 Planetesimal accretion in binaries: latest results 

The stage that has been the most extensively studied in recent years is the one just 
before the final embryos-to-planets phase, i.e., the one leading, through mutual accre- 
tion of kilometer-sized planetesimals, to the embryo themselves. The reason for this 
intense research activity is that this stage is potentially the one that is most aff'ected 
by binary perturbations. Indeed, in the standard planet-formation version, this stage 
proceeds through fast runaway and oligarchic growth that require very low impact ve- 
locities between colliding bodies, typically s maller than their escape velocity, i.e., just 



a few m.s ' for kilometer-sized objects (e.g., Lissauer , 19931) . It is thus very sensitive 
to dynamical perturbations, the crucial parameter sealing the fate of the planetesimal 
population being the distribution of encounter velocities dv among them. As this field 
of research is a (very) fast evolving one, let us here briefly summarize and synthesize 
the main results obtained so far, especially in the past couple of years. 



Early works (e.g. iHeppenheimer 1978 ) revealed that estimating the dv distribu- 



tion is not straightforward. In the complex dynamical environment of a binary, dv 
is indeed no longer directly proportional to the bodies eccentricity e and inclination 
i, because planetesimal orbits are strongly phased and not randomly distributed in 
periastron and ascending nodes. 
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The pioneering study of iMarzari & Scholll (120001) later showed that a fundamen- 
tal mechanism controlling planetesimal dynamical evolution is the coupling between 
secular perturbations forced by the companion and friction with the primordial gas 
that is left in the protoplanetary disc at this sta g e. This cou pling results in a strong 
phasing of planetesimal orbits. iThebault et al. I (l2004 I2OO 6) have shown that this 
phasing is size-dependent, so that dv are small between equal-sized objects but can 
reach very high values for bodies of different sizes. For most "reasonable" size dis- 
tributions within th e planetes i mal p opulation, the differential phasing effect is the 



dominant one (Theb ault et al.L 120061) . This can lead to an accretion-hostile environ- 



ment in vast regions of the circumprimary disc, which can stretch much closer to the 
primary than the radial limit for orbital stability or embryo accretion. Furthermore, 
even in most of the regions where planetesimal accretion is possible, it cannot pro- 
ceed in the same way as around a single star, because the dv increase is still enough 
to strongly slow down and impede the runaway growth mode. A worrying result was 
that, for the emblematic 7 Cephei case, the location at which the planet is observed 
is probably too perturbed to allow for the planetesimal accretion stage to proceed 



(IPaardekooper et al.Ll2008i) . Similarly troublesome results were later obtained for po- 



tential planets in the habit able zone (HZ) of both stars of arguably the most famous 



binary system: a Centauri (IThebault et al.Ll2008Ll2009h . 



These pessimistic studies were all considering a simplified static and axisymetric 
gas disc. However, the simulations of Paardekooper et al. (2008) showed that the situ- 
ation gets even worse when considering a dynamically evolving gas disc: gas stream- 
lines and planetesimals follow very different orbits, which in most cases increases gas 
friction and thus the accretion-hostile effect of differential^phasing of planetesimals. 
On a more positive note, IXie & Zhoul (f2009l) and Ixie et al.l (l2010al) showed that a 
small inclination of a few degrees between the circumprimary gas disc and the binary 
orbital plane could in fact help accretion. This is because planetesimal orbital incli- 
nations are segregated by size, thus favouring low-dv impacts between equal-sized 
bodies over high-iiv impacts between differently-sized objects. As a result, planetes- 
imal accretion could become possible, albe it in a v ery slowed down form, in the HZ 
of a Cen. However, a very recent study bv iFragner et al.l (1201 li) seems to indicate that 
taking into account the effect of the gas disc's gravity could offset this positive effect 
of orbital inclination, leading to high dv dynamical environments. 

This planetesimal accretion issue is thus far from having been solved, and is in 
any case very dependent of the set up: binary orbit, relative masses of the star, location 
in the circumprimary disc, etc. (Thebault et al.. 2006,) . 



1.4 HD 196885 

Given the huge parameter space to explore, most studies of planet formation in bi- 
naries consider only one specific and illustrative example. It is usually either a Cen- 
tauri, because o f its obvi ous interest as our closest neighbour and potential for planet 
detection (Gu edes et al.[ [2008). or y Cephei, because it was until recently the most 
"extreme" exoplanet in a binary, in the sense that it is the one, given its large distance 
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to the primary and proximity to the secondary, for which companion perturbations 
are expected to be the strongest. 

This "privileged" position of j Ceph has however been challenged by a recent 
study bv lChauvin et alj (1201 ih on HD196885 . This system had already been identi- 
fied as a planet-hosting binary (IChauvin et all l2007h . but only the comp anion's pro- 



j ected distance was known, not the binary's orbit. With new observations. iChauvin et al 
(1201 ih were able to constrain this orbit to = 21 .0 ± 0.86 AU and cb- 0.42 ± 0.03, 
and refined the exoplanet orbit to ap = 2.6 + 0.1 a nd ep = 0.48 ± 0.02 (a value 
close to the earher esimate by Correia et al.l 2008 ). This makes HD196885Ab a 
more extreme binary exoplanet than y C ephei Ab. One way to quantify it is by 
looking at the Holman & Wiegertl (1999) criteria for orbital stability, which gives 
a critical radial distance (to the primary) of flcn/ - 3.71 AU for y Ceph and 3.82 AU 
for HD 196885, thus placing the planet at ~ 1.7 AU from the orbital stability limit 
in y Ceph but at only ~ 1.2 AU for HD196885. In fact, this planet's eccentricity, 
ep = 0.48, makes it reach ~ 3.85 AU, i.e. slighly beyond Ucrii, although this does not 
necessaril y mean its orbit is unstable , since Ocrit is not a razor sharp boundary (the 
formula of Holman & Wieger^ ( 1999 ) comes with error bar s delimiting a ''gray" a rea 
around fle„,). However, a preliminary numerical analysis by Chauvin et al. ( 201 11) in- 
dicates that this planet might possibly be on an unstable orbit, unless there is a high 
(unconstrained) inclination between the binary plane and the planet orbit. These sta- 
bihty issues need to be further investigated, but it is clear that this system is by far the 
most "extreme" planet-in-a-binary so far 



1 .5 present work 

We shall not reinvestigate the issue of the planet's long-term stability, nor shall we in- 
vestigate the possible cause for its high eccentricity of 0.48. We shall here focus on the 
key issue of the formation of a planet in such a perturbed environment. HD 196885Ab 
is indeed clearly the planet that poses the strongest challenge to any planet-formation 
model, especially regarding the planetesimal-accretion stage. We numerically investi- 
gate under which conditions this stage might, or might not, proceed in HD196885A's 
circumprimary disc. As with most previous similar studies, we follow the distribution 
of one crucial parameter: impact velocities amongst planetesimals. Our aim is to iden- 
tify the regions where the dynamical environement is too perturbed to allow classical 
core-accretion of a planet. To obtain conservative results, we shall consider the most 
favourable (i.e., accretion friendly) assumptions for our simulation: axisymetric static 
gas disc, no self gravity. We present our numerical model in Sec|2] Our main results 
are presented in Sec|3] In Sec|4] we discuss the robustness of our main conclusions, 
i.e., that most of the circumprimary disc is hostile to accretion, and consider several 
possible solutions to the apparent paradox of having a planet in a accretion-hostile 
region. Conclusions are presented in SecjS] 
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2 Model 



We follow a population of = 2 x 10"* planetesimals, sampling a much larger pop- 
ulation of "real" physical planetesimals, orbiting the primary and dynamically per- 
turbed by the companion. The forces acting on the particles are both stars' gravity 
and friction with the primordial gas disc. The code has a built-in collision search 
algorithm, tracking, at each time step, all mutual encounters between all bodies. In 
order to yield a statistically significant number of encounters despite the limited num- 
ber of particles as compared to a real planetesimal populat ion, we resort to the usual 
method of assigning an inflated radius to each particle (e.g.lTheb ault & Brahic, 1998_ 

Marzari & Schollll2000l;[chamoz et all 12001 1; Imi wick & Chia'na. ,2007: Xie & Zhqu^ 

20081) . For a more detailed descrip tion of our algorithm, see for example T hebault & Brahic 
(Il998h and lThebault et al.1 (l2006l) . 



Gas drag is computed following lWeidenschilling & DavisI (Il985h : 



(1) 



where F is the force per unit mass, Vp_g the velocity of the planetesimal with respect 
to the gas, Vp^g the velocity modulus, and K is the drag parameter given by: 



K 



3pgCd 



(2) 



where pg is the gas density, and ppi and s are the planetesimal density and radius, 
respectively. The coefficient Cd is a dimensionless quantity related to the shape and 
size of the body 0.4 for spherical bodies). The gas disc is assumed to be static and 
axisymetric. For the gas density and gas velocity, we follow iTakeuchi & LinI (120021) 
and assume 



Pg(r,z)=Pgo\ 



AUI 



exp 



2 \ 



Z 



and 



Vg(r, Z) = Vk,mid 



2 r 



g > 



(3) 



(4) 



where hg{r) - /!o(r/AU)**^^'^^ is the scale height of the gas disc and vi,^^^ is the 
Keplerian ve locity in the mid plane. We consider the Minimum Mass of Solar Nebula 
(MMSN, see lHavashiL[l98ll) as a reference disc, where p - -2.75, q = -0.5, pgo = 
1.4 X 10-^g.cm-^ and ho = 4.7 x lO^^ AU. 

The static and axisymetric assumption are taken for computing time reasons but 
is of course a crude simplification of the real behaviour of a circumprimary gas disc, 
which should react to the companion's pert urbations and display pronoun ced eccen- 
tric shapes and azimutal anisotropies (e.g. Artvmowicz & Lubow, '1994*). However, 
preliminary studies with evolving gas discs (Paardekooper et al., 2008) have shown 
that gaseous friction, and in particular the differential phasing effect according to 
planetesimal size, is higher than for a static gas disc. This is because the gas disc gets 
eccentric (Goodchild & Ogilvie, 2006), with an eccentricity eg and a precession rate 
that in most cases strongly departs from that of the planetesimals (see for instance 
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Table 1 Setup for the nominal ran (see main text for parameter definition) 

Binary mass ratio fi = 0.35 

semi-major axis at = 21.0 

eccentiicity ej, = 0.42 

Number of test particles 2 X lO** 
Inflated radius (collision search routine) 5 X 10"^ AU 

Physical radius 1 km< s <10km 

Initial semi-major axis 0.9 AU< a <3.1 AU 

Initial eccentricity < e < 5 X 10"^ 

Initial inclination < (' < 2.5 X 10"^ 

Gas Disc density at lAU po = 1-4 X 10^'g.cm"' 

radial profile Pg('') r"^ 

scale heigh hg = 4.7 X 10-2(r/lAU)' AU 

vertical profile PgU) exp(-z^/2/ij) 



a clear illustration in Figs. 9 and 10 of iPaardekooper et al.i l2008b . This increases the 



relative velocities between planetesimals and gas streamlines, and thus the friction of 
the latter on the former This makes the systems globally more accretion hostile than 
in the axisymetric case (see Sec jl.3b . As a consequence, the axisymetric assumption 
should be regarded as a limiting best-case scenario for planetesimal accretion. 



2.1 setup 

We consider a disc of planetesimals with initial semi-major axis 0.9 < a < 3.1 AU, 
having randomly distributed orbits (longitude of periastron and of ascending node) 
and initial eccentricity < e < 5 x 10"^ and inchnation < / < 2.5 x 10"^. This 
ensures that initial encounter velocities are such as dvi„i, ~ 1 - 2m. s."', approxi- 
mately the escape velocity of a 1km body. This is the velocity distribution expected 
in an unpertur bed population of km-sized planetesimals, in which runaway growth 



can proceed (Lissa uerL 119931) . The planetesimals physical sizes are randomly dis- 
tributed between 1 and 10km. The inflated radius for the collisional search routine is 
5 X 10"^ AU. This size is large enough to yield a statistically significant number of 
impacts, but small enough not to introduce any bias in estimating dv. The set-up for 
our nominal run is summarized in Tab{T] 



2.2 Accretion and Fragmentation Prescription 

Our simulations provide us with the distribution of encounter velocities dv(si,.i2) for all 
impacting planetesimal pairs of sizes si and S2- A key issue is then to interpret these 
velocities and see if they are low enough to allow accretion or are on the contrary too 
high and lead to mass loss or fragmentation of the impactors. Three different regimes 
are possible, defined by two critical velocities Vesc(s\,s2) and Vei-o{si,s2)' 

- dv(si^s2) ^ Vi,sc(si,s2)' "Unperturbed" case. The impact velocity is only marginally 
increased with respect to its initial t/v,„,-, ~ 1 - 2m. s~' value and stays below the 
escape velocity Vesc(s\,s2) for the impacting pair. In this case, "normal", single-star 
like runaway accretion is possible. 
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Vesc(si,s2) ^ dv(si,s2) ^ dvero{si,s2)' "Perturbed accretion" case. The impact velocity 
is increased beyond the escape velocity. This switches off the runaway growth 
mode. However, dv stays at a value small enough to allow some accretion between 
the impacting bodies. 

Vero(si,s2) ^ '^t^'(.!i,.s2): "Erosion". In this case velocities are too high to allow accre- 
tion. Each impacts results in mass loss, i.e., erosion or fragmentation of one or 
both impactors. 



While the value of Vesc(si,s2) is easy to derive, the erosion threshold velocity 
Vero(s\,s2) IS much more difficult to estimate, as it depends on many physical parame- 
ters (particle sizes, compositon, impact velocity and angle, etc..) and because, even 
for the same set of parameters, there exi sts many diverging est i mates of the accre- 



tion/er osion limit (see the discussi on in iThebault & AugereauL l2007h . In previous 



papers (IThebault et al.L 12006. .20081 120091) we considered a complex, and rather cum- 
bersome prescription for Vero(si.s2), trying to connect several possible impact regimes 
(cratering, shattering) and to synthesize sever al available estimates from t he littera- 
ture. In the meantime, a remarkable study by Stewart & Leinhardi (2009) has been 
published, presenting an innovative and simplified criteria for the disruption of plan- 
etesimals. Even if this study's title claims that it is only valid for the "catastrophic" 
disruption of planetesimals, it is in fact applicable to a larger domain, including im- 
pacts usually described as non-catastrophic (where the biggest r emaining fragment is 
more than half the mass of the impactor ). Although the results oP Stewart & Leinhardt 
(12009) have their limitations, we follow iFragner et al.l(l201 1.) and adopt their prescrip- 
tion here because of its simplicity and its self-consistency. Its main parameter is the 
velocity-dependent "reduced" catastrophic disruption specific energy 



Qri 



,9ff/(3-2<P)j^2-3ff 



(5) 



where Sc - (Sj + ^2)'^^ th^ reduced radius of the combined projectile and target 
mass, and a and are material properties. If Qr - Q .5m\m2dv j {m\ + 1112)^ is the 
reduced kinetic energy, then the mass m/^ of the largest remaining fragment is given 
by 



mir 



(mi + 1112) 



1-0.5 



Qr_ 



(6) 



RD 



With the convention that nii > m2, the criteria for accretion is then niir > nii, 
which translates into 



dv < V 



erois 1 ,^2) 



0.5 



(7) 



We follow [stewart & Leinhardt (2009) and consider two limiting cases: v^roi for 
weak rocky aggregates (a - 0.4, - 7,qs = 500, qg - 10^'^, in cgs units) and Vero2 
for strong compact rocks (a = 0.5, <t> -d>, qs -1 x 10^*, qg = 10""^). 
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1.0 1.5 2.0 2.5 3.0 1.0 1.5 2.0 2.5 3.0 

semi-major axis (A.U.) distance to the primary (AU) 



Fig. 1 Gas free case. Snapshots, after t= 4 X 10^ and t= 2 X 10** years, of the eccentiicity distribution 
(left-hand side) and encounter velocity distribution (right-hand side) within the planetesimal population, 
as a function of semi-major axis (for the eccentricity) and radial distance to the primary {dv panels). The 
straight line on the (e,a) panels indicates the value of the forced secular eccentricity ey. The velocity 
distribution is obtained by recording all impacts in the t± 100 years time interval. The inward propagation 
of the high-dv "wave", due to orbital crossing of neighbouring orbits, is clearly visible. 



3 Results 



3.1 Gas Free Case 



In order to clearly identify the different mechanisms at play, we first present a fiducial 
and pedagogical run with no gas. Note that planetesimal physical sizes are irrelevant 
for this purely gravitational case. Figs{TJ and c show the classical build up of large 
secular eccentricity oscillations around the forced eccentricity e j. These oscillations 
are due to the fact that, at each radial distance from the primary, eccentricities vary 
with time following a sinusoidal function 



e{a, f) -Icf 



5 a 

sin(ut/2) \- — 



2atl - el 



sin(ut/2) 



(8) 



whose frequency u vary with semi-major axis (e.g. iThebault et al.L l2006h . As has 
been pointed out in previous studies, such oscillations do not immediatly lead to high 
impact velocities, because neighbouring orbits are strongly phased. However, because 
the a dependency of particle eccentricities increases with time (since u depends on a), 
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the e(fl)-oscillations get narrower with time and orbits within one oscillation "wave" 
eventually c ross, at which point very high dv are suddenly reached (see discussion in 



Thebault et al- .2006) . The radial location at which orbits cross m oves inward with 



time following the empirical law derived bv lThebault et al. (l2006h 

/ .n1.07 



^-^'■''^^[ul] (loXu) (lov) 

Note that the location and timing of the orbital crossing does only weakly depend 
on the initial conditions for planetesimal orbits. We have chosen here the simplest 
and probably les s unlikely case of ini tial circular orbits (for more on this issue, see 
the discussion in Thebault et al. but taking initial eccentric orbits will roughly 



lead to the same behaviour for across, the only change being an additional free com- 
ponent to the encounter velocities. 

From FigllJ), we see that across, i-C- the location beyond which there is an abrupt 
increase of impact velocities, reaches 2.6 AU (the present location of the planet) in 
less than 4 x 10^ years. At this point, encounter velocities increase by more than a 
factor ~ 20. However, even the short timespan before orbital crossing is not fully 
calm, because of sporadic high-iiv impacts (100 to 200m. s"') with particles in the 
nearby mean motion resonances clearly seen on FigH^. In this gas-free case, the 
r ~ 2.6 AU region is thus very hostile to low-dv planetesimal accretion. 

The situation is less desperate closer to the primary. The regions inside 1.5 AU is 
for instance protected from orbital crossing, and high-t/v, for more than 2 x 10"* years 
(FigHJl). This should in principle leave enough time for runaway accretion to produce 
embryos from kilometre-sized planetesimals. However, we shall see that gas drag 
completely obliterates these optimistic conclusions. 



3.2 Nominal Run with Gas 

Gas drag radically changes the dynamical evolution of the system. Its main effect 
is to phase both planetesimal eccentricities (Fig|2^ and c) and longitude of perias- 
tron 0) (Fig|2) and d) according to their size s. As shown in Fig|2};, the initially 
large eccentricity oscillations are progressively damped. The system will eventually 
tend towards a steady state, where all eccentricities reach an equilibrium value e^a^ ,) 
depending on semi-major axis and size. In the innermost regions with higher gas den- 
sities, this steady state is reached relatively quickly for all particles in our 1-10 km 
size range. Beyond ~ 1.4AU, however, the steady state has not been reached, at 
f = 2 X 10'* years, for the biggest planetesimals. And at the location of the planet, 
2.6 AU, even the smallest 1km objects still have residual eccentricity variations at 
the end of the run (Fig|2};). 

These dynamical behaviours have clear consequences on impact velocities, which 
reach high values everywhere in the disc (Figj3). These high dv are reached very 
quickly, a few 1000 years, and this for two different reasons depending on location 
within the disc: 
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Fig. 2 Nominal case with gas friction from a IxMMSN disc: eccentricity (left panels) and longitude of 
periastron (right panels) as a function of semi-major axis (the binary's longitude of peiiastron is 0), at 
f = 4 X 10^ and r = 2 X lO" years, for a population of planetesimals in the 1 < .v < 10 km size range. 




Fig. 3 Same nominal gas-friction run as in Fig|2] impact velocity distribution as a function of radial 
distance to the primary. The colours indicate the expected collision-outcome regime for each impact: "un- 
perturbed" runaway accretion (green), perturbed accretion (blue), erosion {red). The yellow area is for 
impacts where v^rol <dv< Vero2, for which the collisional outcome is uncertain (see text for details). 
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Fig. 4 Nominal gas drag case. Time evolution of the impact velocity distribution at 2 different locations 
in the disc. All impacts are recorded in a narrow ring of width 0.02 AU around the central location. The 
colour scale is the same as in Fig[3] 



- In the inner disc, shortwards of ~ 2-2.5 AU, where gas densities are high, the 
differential phasing induced by gas drag, and the velocity increase that comes 
with it, is felt very early on, long before orbits reach a steady state with fixed e and 
(1). This is illustrated in Fig|3^ showing that, at f = 4 x lO^years, when no region 
of the disc has reached a steady state yet, dv already reach values > 300 m.s"' in 
the whole < 2.5 AU region. 

- In the r > 2.5 AU region, the gas drag-induced dv increase is weaker, but veloc- 
ities nevertheless reach even higher values, ~ lOOOm.s"', because of the purely 
dynamical orbital crossing effect identified in the gas-free runs. This is clearly 
illustrated by the dv "jump" at ~ 2.5 AU in Fig[3^. 

. The velocity distribution at f = 2 x 10"* years is remarkably similar to the one at 
4 X lO-' years. The only difference is in the outermost regions, where no dv jump 
is longer visible, because gas drag is now the dominant dv-inducing mechanism in 
the whole system Q (even if a steady state has not been reached in these outer re- 
gions). To illustrate these behaviours more clearly, we display in Fig|4]the temporal 
evolution of dv in two opposite regions of the disc. At 1 AU, differential phasing 
induces a high-t/v regime after ~ 2 x lO^years, while a steady state is reached af- 
ter ~ 5 x lO^years. At 2.6 AU (the planet location), three succesive phases can be 
distingued: a first phase, starting after only a few 100 years, of moderate-to-high-c/v 
induced by sporadic impacts with objects in neighbouring resonances, followed by 
a second stage, at f = 4 x lO-'years, when orbital crossing occurs and increases dv 
to even higher values, and finally a third stage, starting around 1.5 x lO^years, when 
gas drag phasing progressively takes over as the dominant c/v-inducing mechanism. 
Despite these differences, however, the important result is that a high-t/v regime, re- 
gardless of its different causes depending on location in the disc, is reached after only 
a few 1 000 years in both the inner and outer regions of the disc . 

To estimate the global consequences of this high-t/v regime on the accretional 
evolution of a "real" population of planetesimals, one has to consider a parameter 

' This f ~ 2 X lO'* years time is thus the characteristic timescale for gas drag to dominate the whole 
system's dynamics 
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that has been ignored so far: the planetesimal size distribution. In our simulations, a 
flat distribution between 1 and lOkm has been considered for the sake of simplicity, 
but real distributions should be more complex. This issue of the initial planetesimal 
size distribution is a difficult one. In models of planetesimal a ccretion in our solar 



system, a singl e "initial" planetesimal size is usually assumed dMakino et al.L 11998 



Kokubo & IdaL boOO), but this assumption is taken for the sake of simplicity and 



there should be some size dispersion in any realistic initial planetesimal population. 
The exact profile of this initial distribution is difficult to constrain. Firstly because it 
depends on the way planetesimals are formed from smaller grains and pebbles, a pro- 
cess which is still far from b eing fully understood, even if significant p rogresses have 
been made in recent years (I johansen et al. , 2007t Cuzzi et al. , 2008i) . Furthermore, 
the very concept of an "initial" size distribution can be questioned, as there might be 
a wid e spread in the times at whi ch km-sized objects appear in a given region of the 
disc (iChambersl I2OIOI: IXie et al.L 1201 Ob . These iss ues go well bey ond the scope of 
the present paper, a nd we shall consider, following iThebault et al. I {2OO8, 2009) and 
Xie & Zhoul (I2OO9I) . a Maxwellian distribution centered on i = 5 km. Such a rela- 



tively peaked disitribution is, in line with our conservative approach, a priori more 
accretion-friendly since it minimizes the rate of encounters between differently-sized 
objects. It also agrees with most planetesimal-formation scenarios' conlusion that 
there should be a privileged size for initial planetesimals. In practice, we weight each 
S1-S2 impact obtained in our run with a flat size distribution by a factor /(.5,,j2) ac- 
counting for the Maxwellian distribution 0. 

The accretion/erosion behaviour of the whole disc is displayed, at f = lO'^years, 
in Fig|5] As can be clearly seen, the rate of impacts leading to mass erosion or frag- 
mentation is > 80% everywhere, except around 1.6 AU where it is ~ 65%. When 
discarding the impacts with "uncertain" (yellow) outcome, the level of accreting im- 
pacts is less than 10% everywhere. Moreover, among these accreting impacts, most 
of them are in the "perturbed" mode (blue), with almost no impact allowing runaway 
growth (green). We have tried different size-distributions and always found the same 
global accretion-hostile trend, except for extremely peaked, and probably unrealistic, 
distributions. Only with much larger planetesimals can this negative trend be reversed 
(see SeclX4li. 



3.3 Inclined Binary 



Xie & Zhou (l2009h have shown that a small inclination between the binary and the 



circumprimary disc can help accretion by segregating particle inclinations accord- 
ing to sizes, thus favouring impacts between equal-sized bodies. We explored this 
possibility for the present HD 196885 case, assuming a small inclination of 2" for the 
binary. As shown in Fig|6] an improvement is obtained compared to the coplanar case 
(Fig©. Accreting impacts now make up between 5 and 20% of all impacts in most 
of the disc. However, the system is still globally hostile to accretion everywhere (be- 
tween 60 and 80% of "red" impacts), especially at the location of the planet (2.6 AU) 



- The reason why we do not run a simulation with a Maxwellian distribution to start with is because, 
with a flat 1-lOkm distribution, we can get a good statistics on all possible impacting sizes si — S2 
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Fig. 5 Nominal gas drag case. Relative importances, att= 10 years, of the 4 possible collision outcomes 
as a function of radial distance to the primary (the colour scale is the same as in Figs|3]and in Fig|4). A 
Maxwellian size distribution centered on 5 km is assumed. 




radial distance (AU) 

Fig. 6 Same as Fig|5] but with a binary inclined by 2" with respect to the circumprimary disc. 

where the fraction of eroding impacts is in excess of 90%. We have tried several dif- 
ferent values for the binary inclination in the 1" < / < 10" range, and always end up 
with results very similar to those displayed in Figi6]i.e., a system that is globally very 
hostile to accretion. 

Strengthening our conclusions is also the fact that the vertical segregation accord- 
ing to particle sizes obtained bv .Xie & Zhou (2009) is probably unrealistically high. 
This is because it has been obtained for an axisymetric gas disc, whereas a real gas 
disc would get eccentric and tend to diminish the size-sorting effect. Our results for 
an inclined binary, also obtained with an axisymetric gas disc, should thus be con- 
sidered as a best-case scenario regarding planetesimal accretion; a best-case scenario 
for which no accretion is possible in the whole > 0.9 AU region. 



3.4 Large planetesimals 



As suggested bv [ Thebault et alJ (120081) and, using slighly different assumptions, by 
Beauge et al. I (1201 oh . another possible solution to the accretion-hostile-velocities dilemma 
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Fig. 7 Large planetesimals (10 < s < 100 km) case. Encounter velocity distribution att= 10 years (left) 
and respective balance between collision outcomes (right). The colour scale is the same as in Fig|5] 



is to start from larger planetesimals. We explore this hypothesis by considering a pop- 
ulation of initial planetesimals in the 10-100 km range. Taking larger objects reduces 
the impact of gas drag and its subsequent effect on increasing dv among planetesi- 
mals. This is what is observed in the inner regions of the disc in Fig|7h, where impact 
velocities are lower than for the 1-lOkm population (Figj3^). This situation is more 
favourable for accretion, a tendency that is amplified by the fact that, for a given dv, 
bigger planetesimals are more resistant to impacts than smaller ones (in the gravity 
regime valid for objects in the > 1-lOkm range). As a result, at f ~ 10"* years, the 
whole region shortwards of ~ 2 AU is accretion-friendly for 100km planetesimals 
(FigGh)- In the outer regions, however, the situation is radically different: the values 
of dv are much higher because of the secular orbital crossing effect. In these outer 
gas-poor regions, large planetesimals behave almost as test particles in the gas-free 
case. In regions of orbital crossing, impact velocities are high enough to lead to ero- 
sion even for 100 km objects (Fig|7j)). As in the gas-free case, the high-c/v orbital 
crossing "front" reaches 2 AU in ~ 10"* years, but the 2.6 AU location is reached in 
less than 4 x 10^ years. 



4 Discussion 



The previous results all tend to indicate that the HD196885 circumprimary disc is 
strongly hostile to the accretion of kilometre-sized bodies. This is especially true of 
the region at 2.6 AU from the primary, i.e., around the current location of the detected 
planet. Planetesimals are basically caught between a rock and a hard place: in the 
inner regions secular perturbations combined to gas drag induce high-iiv between 
all non-equaly-sized objects, while in the outer regions high-t/v are due to secular 
perturbations alone, which make orbits cross within a few thousand years. 



16 



P. Thebault 



4. 1 Limitations, robustness of our conlusions 



As underlined earlier, our numerical exploration relies on several simplifications. It 
has also many free parameters that cannot be all thoroughly explored. How robust 
are our conclusions, especially regarding the 2.6 AU region, with respect to these 
limitations? 



Timescale. We first note that timescale is not a critical issue: high-t/v are reached 
almost everywhere after only a few lO^years. The only change with time is that gas 
drag progressively takes over as the dominant c/v-inducing process even in the outer 
regions. 



Binary inclination. As shown by IXie & Zhou (l2009h . a small inclination between 
circumprimary disc and binary orbital plane acts in favour of planetesimal accre- 
tion. This point is all the more appealing because this inclination is in most cases 
an unconstrained parameter, even for close binari es, for which the assumed "copla- 
narity" is not constrained to less than ~ 10° dHalel 1 19941) . In the present case, 
however, even if the situation improves significantly for ib of a few degrees, it is far 
from enough to reverse the general accretion-hostile trend of the system (Fig|6ll. In 
addition, as pointed out in Sec j3.3l the vertical size-sorting of planetesimals, and its 
effect in favour of accretion, is probably overestimated in a circular gas-disc case. As 
a consequence, we believe our conclusions of an accretion-hostile disc with a slighly 
inclined binary to be relatively robust. 

Of course, even higher values of ih c ould be possible, as m ight be suggested by the 
preliminary long-term stability study of lChauvui et al.l(l201 U . These high inclination 
cases cannot be explored with the present model, and will be the purpose of a forth- 
coming general study devoted to this issue (Xie et al., submitted), but it seems very 
unlikely that impact velocities should be reduced in such high-//, systems, especially 
when the Kozai regime sets in. 



Gas disc profile. We have run several additional simulations (not shown here), ex- 
ploring different gas disc profiles and densities. In almost all cases, results are roughly 
comparable to the nominal case: high, and accretion-inhibiting impact velocities in 
the whole r > 0.9 AU region. Only for very tenuous discs do we get an accretion- 
friendly inner region for more than lO^^years, which is basically the gas-free result 
displayed in FiglT] It coud be argued that this gas-free stage is the one towards which 
the real circumprimary disc is naturally evolving, since primordial gas is expected to 
be removed after a few 10^ years in protoplanetary discs, or even less than that for 
close binaries dCieza et al .I I2OO9I) . In fact, iXie & Zhoul (l2008h have shown that some 
orbital rephasing could occur during the gas removal phase, possibly rending the disc 
accretion-friendly again. However, as shown by Thebault et al. (2008), it is unlikely 
that planetesimals could survive the long accretion-hostile period before gas removal 
without being grounded to dust and removed by gas friction. The only solution would 
be that planetesimals form late in the disc's history, when few or not gas is left, but 
this hypothesis conflicts with all planetesimal formation models. In any case, let us 
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Stress that even in this unlikely scenario, the region around 2.6 AU would still be 
hostile to planetesimal accretion, due to secular effects alone. 

There is however another important gas-disc-related issue to consider here, i.e., 
that all our runs follow the same crude simplification, inherent to our approach, of a 
static axisymetric gas disc. Nevertheless, as discussed in Sec|2] we expect planetes- 
imals imbedded in evolving g as discs to have impact ve locities that are even higher 



than in the fiducial static case (iPaardekooper et al. So here again, our dv esti 



mates do give a conservative lower estimate. 

Gas self gravity. Another simplification of our model, as well as of most previous 
studies of planetesimals in binaries, is the neglect of the gas disc's grav ity. To our 



knowl edge, the only published studie s taking into account disc gravity are Kiev & Nelson 



(l2007h and the very recent work bv iFragner et al.l (1201 11) . These pioneering studies 



present results with a limited number of particles (which is the price to pay for in- 
cluding the disc's gravity) that does not allow accurate dv estimates, and for which it 
is difficult to untangle the effect of gravity from that of gas drag. However, it appears 
clearly that the global qualitative effect of disc gravity is to further increase impact 
velocities, by adding an additional jitter to the eccentricity and periastron evolution of 
planetesimal orbits |j. We do thus expect our gas drag-only simulations to here again 
give a lower limit for "real" impact velocities. 



4.2 A planet in an accretion-hostile environment? 



That most of the circumprimary disc is too excited to allow planetesimal accretion 
seems to be a relatively robust result for the HD 196885 system. And yet there is 
a planet well inside this accretion-hostile region. There are basically four potential 
solutions to this paradox: either 1) the planet could form in situ by being able to 
bypass the mutual-planetesimal-accretion phase, or 2) it was formed elsewhere and 
was later injected at its present location, or 3) the binary had a different orbit during 
its early history, or 4) the planet did not form by core-accretion but by direct disc 
instability. 



4.2.7 Bypassing the kilometre-sized planetesimals accretion phase': 



The most obvious potential way to bypass this stage is if planetesimals were formed 
big, i.e., not in the kilometre-sized range but rather in the 50-100 km one. Interest- 
ingly enough, this big-initial-planetesimals scena rio seems to be the on e fa voured by 
the mo st recent planetesimal-for mation models oflJohansen et alj (120071) and lCuzzi et al 
(120081) . Moreover, according to Morbi delli et al. (l2009l). there seems to be observa- 
tional evidence for an initial population of large, > 100 km, bodies in the asteroid 
belt, even if this conclusio n h as been recently q uestionned, for different reasons, by 
Minton & Malhotral (l2010l) or lXie et al.l ( 1201 017) For the present problem, however, 
there are some issues with this large-planetesimals hypothesis. The first one is that, to 



' This is also the preliminary conclusion of iMarzari et al.l ( I2008I) . who also investigated disc gravity, 
but in the different context of a circumbinary disc. 
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our knowledge, there has been no study of how these big-planetesimal formation sce- 
narios, for which several problems remain to be solved even for a normal single-star 
environment, could proceed in the highly perturbed environment of a close binary. 
The second, and more problematic issue for the specific HD 196885 case is that, even 
// initial planetesimals are ~ 100 km big, the region around 2.6 AU is still hostile to 
accretion (see the "large planetesimals" run in FigjT). We ran an additional test sim- 
ulations and found that the 2.6 AU region becomes accretion-friendly only for bodies 
with sizes > 250 km. It is far from being assured that planetesimals can be born this 
big, especially in a highly perturbed close binary. 

Another possibility to overcome the mutual-planetesimal-accretion hinder, even 
in the case of small planetesimals, is th e so-called "sno wball" growth mode first en- 
countered in the simulations of Paard ekooper & Lein hardt (2010) and investigated 
in more detail bv lXie et all (l2010bh . In this scenario, planetesimals grow preferen- 
tially by sweeping up of small dust particles, provided that the local mass density 
of solids contained in dust, ps(dust), exceeds that contained in planetesimals, ps(pian.)- 
The snowball growth mode is especially appealing for dynamically excited systems 
such as close binaries because the accretion of dust onto planetesimals should toler- 
ate much higher velocities than the mutual accretion of the planetesimals themselves. 
Exactly how much higher is not clear yet, as there exists to our knowledge no pub- 
lished study of the velocity dependence of the accretion efficiency of dust on large 
targets, at least for the dv > 500m. s ' regime encountered hereQ One additional 
issue is of course whether or not the ps(dust) ^ Psipian.) criteria is met in real sys- 
tems. As shown bv lXie et all (l2010bh . this condition could be fulfilled if there is a 
significant spread in the times at whi ch planetesim als do appear in the system, as 
would for instance be expected for the Cuzzi et aP (l2008l) scenario. But even if that 
is not the case, the ps{dust) ^ Psipian.) condition could be met later on in a perturbed 
system's evolution, when the quantity of dust produced by destructive planetesimal- 
on-planete simal impacts exceeds the mass left in the remaining unshattered plan- 
etesimals ( Paardekooper & Leinhardd I2OIOI) . Both of these cases have so far been 
studied in preliminary works using simplified prescri ptions : analyt ical expression for 
the growth of initially isolated planetesimals in Ixie et al. I (1201 Obi) , and 2-D simula- 
tions for the [Paardekooper & Leinhardtl(l2010 l) studies of planetesimal re-accretion of 
impact-produced dust. They both will be quantitatively re-investigated in forthcom- 
ing studies (Paardekooper et al., 201 1 and Xie et al., 2011, both in preparation). 



4.2.2 Embryo migration, Planet-Planet Scattering? 

Apart from the strongly revised versions of the planet-formation scenario presented 
above, there exist other potential solutions to the inhibition of the planetesimal-to- 
embryos stage. The first one is that embryos form in accretion -friendly regions closer 
to the primary and later migrate outward to the planet's present position, where they 
can continue to grow because the final embryo-to-planet stage is much less affected 
by binary perturbations (see Sec{T}. This scenario has been quantitavely investigated 

^ The "high veloc ities" dust impacts considered in laboratory experiments such as those of 
iTeiser &Wurm|j200^ do not exceed 50-lOOm.s 
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bv iPavne et al. I (l2009h . who showed that a fraction of the embryos formed in the 
inner regions can indeed later move out. The relative amplitude 5a I a of this outward 
migration can reach 0.3 to 0.8. This is however not enough for the present case, 
because it means that the innermost possible origin for an embryo having moved to 
2.6 AU is ~ 1.4 AU, i.e., a region that is still highly hostile to planetesimal accretion 
(see Fig|5j- 

On a related note, one could imagine that the planet fully formed in the inner, 
accretion friendly regions and was later ejected by gravitational interactions with a 
second, yet undetected, planet. There is however a major problem with this scenario, 
which is that the accretion-friendly region around HD196885A is very narrow. By 
running an additional simulation focusing on the inner r < 0.9 AU disc, we indeed 
find that the limit between the eroding (majority of "red" impacts) and accreting (ma- 
jority of "green" + "blue" impacts) regimes is located at around 0.4 AU from the 
primary. This would mean that 2 giant planets would have to form within 0.4 AU 
from HD 196885 A. Such an hypothesis seems to be ruled out by all planet-formation 
scenarios. Another counter argument is that a second planet closer to the primary 
would have been detected in radial velocity measurements (Chauvin, personal com- 
muinication). 



4.2.3 Wider initial binary? 

Another possibility is that, during these early stages of its existence, the binary's or- 
bit was different from what it is today. This might happen because most stars are 
expected to be born in clusters, which are initially compact, thus allowing frequen t 



interactions between neighbouring stars. The simulations of Malmb erg et al.l (120071) 



have shown that, for the typical cluster they considered, binaries with moderate-to- 
large separations do suffer early encounters that have on average shrinked their or- 
bit. Interestingly, HD 196885 is just at the limit present-day separation, ~ 20 AU, 
above which these orb it-shrinking effects are found to be significant (see Fig.4 of 
Malmberg et al. . l2007h . There is thus a non-negligible chance that its initial orbit. 



during the planet formation phase, was wider than what it is today. How much wider 
would it need to be in order to allow planetesimal accretion at 2.6 AU? We ran a se- 
ries of test simulations with increasing values of the binary separation and found that 
the r ~ 2.6 AU region becomes accretion-friendly for ab > 45 AU (for the same value 
of the eccentricity e/, = 0.42). This means that the binary's orbit would have had to 
shrink by at least 24 AU during its early history. Is such a value realistic? Unfortu- 
nately, it is impossible to judge from the 'Malm berg et al. I (12007 ') simulations, since 



no information about the amplitude of orbital compaction was presented in this study 
(this important issue should clearly be investigated in future studies). However, such 
a large, ~ 20-25 AU, change in semi-major axis does a priori appear unlikely for a 
binary that is just at the limit separation for which orbital change become significant. 



^ the second undetected planet would also need to be massive in order to be able to perturb HD196885b 
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4.2.4 Formation by disc instability? 



If none of the aforementioned solutions works, then a more radical alternative might 
be considered, i.e., to for go the core-accr etion scenario altogether This is the conclu- 
sion recently reached by Duchenel (12010 ). who argued that the shorter disc lifetime 
in tight binaries makes it difficult to form giant gaseous planets. This study advocates 
a violent formation process, by direct disc fragmentation, for planets in a/, < 100 AU 
binaries, a hypothesis that seems to be supported by the fact that exoplanets within 
ah < 100 AU binaries are significantly more massive than those within wide binaries 
or single stars. Such a scenario is usu ally invoked to ex plain the f ormation of planets 
at large, > 100 AU semi-major axis ( Bolev . 2009h . but Duchenel (12010) argues that, 
in the specific context of tight binaries, it could also work for planets much closer to 
their star His main argument is that circumprimary discs in close binaries could be 
more compact and denser than discs around single stars 0, thus potentially favouring 
gravitational instability, and also because pert urbations of the close stellar companion 
could give an additional trigger to instability (lBossll2006h . 

However, this alternative explanation should be considered with some caution. 
Several studies have indeed shown that the instability scenario does also encounter 
major difficulties in the context of close binaries, a nd that no circ um primary disc 
ets dense enough to be unstable. As an example, iNelsonI (l2000h or iMaver et al. 



2005) conclude that instability is severely hampered by the presence of a close 
(< 50 - 60 AU) companion. This issue is thus far from being settled yet, and fur- 
ther, more detailed investigations are clearly needed to assess if disc-fragmentation 
can be considered as a viable alternative formation channel in close binaries. 



5 Summary and Conclusions 

The planet of the HD196885 system is to this day the one that has the most perturbed 
orbit amongst all the ~ 100 exoplanets detected within binaries. Leaving aside the 
issue of the long term stability of its orbit, the question of how this planet could form 
under such extreme conditions is a critical one, that might have implications on our 
understanding of the planet-formation process in general. 

We have investigated one specific leg of the planet-formation process, the inter- 
mediate stage leading from kilometre-sized planetesimals to protoplanetary embryos, 
as it is the one that is probably the most aff'ected by perturbations in close binaries. 
We numerically follow the evolution of one crucial parameter controlling the accret- 
ing fate of the system: the distribution of impact velocities dv amongst planetesimals 
orbiting HD 196885 A. Due to stringent computing time constraints, we adopt a deter- 
ministic model based on several simplifications (static gas disc, no disc gravity, etc.), 
all of which concurring to give a conservative lower limit for the dv distribution. 

We find that, for almost the whole explored parameter space (planetesimal sizes, 
gas disc profile, inclination between the binary and the circumprimary plane, etc.), 

' lOuchene (2010) argues that the lower (sub)minimetre fluxes measured for discs in < 100 AU binaries 
do not prevent these discs from being as massive (and thus dense r) as sing l e-star ones, because such 
truncated discs get optically thick to their own emission (see Fig. 2 o f lDuchend fSoToh ). 



Planet formation in HD 196885 



21 



impact velocities reach values that lead to eroding impacts between planetesimals. 
Most of the circumprimary disc is thus strongly hostile to accretion, especially the 
region at 2.6 AU from the primary corresponding to the current location of the planet. 

We considered the possibility that the planet was formed in an accretion-friendly 
region much closer to the primary and later moved outward, either during the embryo 
accretion phase or by mutual perturbations with another (undetected) giant planet. 
However, we rule out this hypothesis because the inner accretion-safe region is much 
too narrow and close to the primary, r < 0.4 AU, to allow the formation of two giant 
planets. 

In the highly perturbed r ~ 2.6 AU region, the only way for the HD196885Ab 
planet to form through the mutual accretion of planetesimals is either 1) if these 
planetesimals were initially larger than ~ 250 km, or 2) if the binary had an initial 
separation ai, > 45 AU and shrinked by at least a factor 2 during its early history. 
Although large initial planetesimals as well as early orbital compacting of binaries 
both make sense in view of recent planetesimal-formation and stellar cluster theories, 
the values we find for these 2 parameters do appear rather extreme. An alternative 
solution is to suppose that planetesimal grow through the so-called "snowball" growth 
mode (Xie et al., 2010b) by progressively sweeping up small dust particles. However, 
the efficiency of this alternative growth channel has not been quantitatively estimated 
yet, especially in the specific context of close binaries. 

These results strengthen HD196885Ab's status as the most "extreme" planet-in- 
a-binary known to date. What sets it appart from other cases, like the y Cephei planet 
or a putative habitable planet in the a Cen system, for which the canonical core- 
accretion scenario also encounters serious problems, is that it is the only one for 
which these problems could be serious enough to become insurmontable. Should this 
be the case, then alternative planet-formation scenarios, such as direct collapse due 
to disc instability, migh t be considered for this planet, and possibly all other giant 
planets in close binaries ( Duchene , 2010l) . Unfortunately, the disc instability scenario 



does also encounter severe difficulties in the context of close binaries, and it is too 
early to know it can be a viable alternative formation channel. 
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